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Abstract 

The probability density of low-lying eigenvectors of the hermitian Wilson-Dirac op- 
erator H(k) = 75-Dw( K ) is examined. Comparisons in position and size between 
eigenvectors, topological charge and action density are made. We do this for stan- 
dard Monte-Carlo generated SU(3) background fields and for single instanton back- 
ground fields. Both hot and cooled SU(3) background fields are considered. An 
instanton model is fitted to eigenmodes and topological charge density and the 
sizes and positions of these are compared. 
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1 Introduction 



It is known [1,2,3,4,5] that the zero modes of lattice Dirac operators are 
strongly localised. We show that not only zero modes but all low-lying modes 
of the hermitian Wilson-Dirac operator if (mo) = 7siJw( — m o) are strongly 
localised for bare quark masses < mo < 2. Furthermore, we show that this 
localization is strongly correlated to topological objects, including instantons. 
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A similar examination has been done for the overlap formalism [6]. We are 
interested in examining the extent to which similar properties are already 
manifested at the kernel level of the overlap formalism. 



Our Wilson-Dirac operator is defined in the standard way by 



+(1 + ^)Ul(x - ^(x - fj,)}, 



(1) 



where the hopping parameter k is related to the bare mass by 



k = (8 — 2mo) 



-i 



(2) 



Therefore the above mentioned bare mass range is equivalent to 0.125 < k < 
0.25 at tree level. We examine eigenmodes in and just outside the region m c < 
m < 2, where m c is the "critical-mass" which is at tree-level but for non- 
trivial gauge fields shifts away from 0. This is the range of the mass parameter 
mo needed for use in the overlap formalism. Modes of H{mo) crossing zero in 
this region are accompanied by the abrupt appearance of exact zero modes of 
the overlap-Dirac operator at the corresponding m value [7]. 

We solve the eigenvalue problem D^ip{x) = \ip(x) for the first four low- 
lying eigenvalues. We are interested in calculating eigenmodes in the physical 
region of the overlap formalism with some further points just at the edge of 
this region. This corresponds to a K-region starting slightly below the critical 
value of k = k c and extending slightly beyond the point where doublers appear 
in the continuum limit of the overlap formalism. We consider 0.115 < k < 0.26 
at tree level. Since the critical k- value shifts from its free field value of 0.125 to 
some higher value for non-zero gauge coupling, we have to adjust our K-range 
accordingly. We look for any change in behaviour of the eigenmodes at the 
border of the region of physical interest. 

Eigenmodes are found by an accelerated conjugate gradient routine [8] which 
is further improved by using dynamic state renormalisation. The major ad- 
vantage of a conjugate gradient algorithm besides its almost perfect parallel 
structure is that it yields not only eigenvalues with appropriate degeneracies, 
but eigenvectors, as well. For selected K-values we also calculate up to 20 
eigenmodes. In the following, the phrase "low-lying eigenmodes" should be 
understood to mean eigenmodes corresponding to the low-lying eigenvalues. 

In order to examine localisations of calculated eigenmodes we plot the proba- 
bility density p(x) =|| ip(x) || 2 for three dimensional cuts through the lattice. 
For comparison we plot the action density as well as the topological charge 
density for the appropriate configuration in the same way. Our calculations are 
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made on 8 3 x 16, 16 4 and 16 3 x 32 lattices with anti-periodic fermion boundary 
conditions in the time direction. On the 8 3 x 16 lattice we test the correlation 
of low-lying eigenmodes with a single instanton configuration and standard 
Monte-Carlo generated SU(3) background fields. For the latter background 
fields we consider hot, 5-sweep and 12-sweep cooled configurations. Calcula- 
tions on the larger lattices are to verify the conclusions we obtain from the 
smaller lattice. In order to quantify our results further we fit an instanton 
model to the obtained data. We distinguish between the model for action or 
charge densities and the model for the zero- mode density [9], 



2 p 2 

p{ x )zcto = c • — • — — — —3, (4) 

7r ((x - x y + p 2 ) 



where x is the distance from the instanton peak at xq to the calculated den- 
sities. The normalisation factor c allows us to fit to the instanton shape and 
prevents the fit from being dominated by the maximal value of the fitted object 
which is affected by periodic images due to the finite volume of the 4-torus. 
Both models are continuum results derived from the standard 't Hooft ansatz. 
Eq. (3) is the action density and is used to fit action and charge densities. 
Eq. (4) is the density of the fermion field in the zero-mode and is used to fit 
eigenmode densities. This allows us to compare localisation sizes and positions 
quantitatively. 

This paper is organised as follows: After a discussion in Sec. 2 of lattice tech- 
niques and general properties, we describe results for a single instanton con- 
figuration on a 16 4 lattice in Sec. 3. In Sec. 4 we investigate localisations of 
eigenmodes of the hermitian Wilson-Dirac operator and the correlations those 
eigenmodes have with topology on standard Monte-Carlo generated SU(3) 
configurations on a 16 3 x 32 lattice. We do this for both hot and 12-sweep- 
cooled configurations. We quantify those localisations and their correlations 
to topology and compare sizes of the eigenmode localization with sizes of the 
corresponding topological object in Sec. 5. Conclusions are presented in Sec. 6. 



2 Lattice techniques 



We carry out our studies on three different lattices. Single instanton config- 
urations are created on an 8 3 x 16 lattice with p = 1.0 and on a 16 4 lat- 
tice with p = 2.0. The latter configuration is cooled for 2 sweeps using the 
standard Wilson action to minimize boundary effects. We consider three stan- 
dard Monte-Carlo SU(3) configurations generated on an 8 3 x 16 lattice with 
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/3 = 4.38, which corresponds to a lattice spacing of 0.165(2) fm. The second 
lattice size used is 16 3 x 32 with f3 = 4.60 which corresponds to a lattice 
spacing of 0.125(2) fm. These configurations are generated with a plaquette 
plus rectangle improved action with mean-field improved coefficients. For the 
smaller exploratory lattice the standard Wilson action and standard plaquette 
topological charge are used to minimize boundary effects. For the bigger lat- 
tice, a three-loop improved action and three-loop improved topological charge 
operator [10,11] are used. 

The topological charge density is given by 

g 2 

l( x ) = ^r^^^(F„u(x)F pa (x)). (5) 



For the standard plaquette topological charge we use the traceless definition 
of F^v extracted from the consideration of 1 x 1 plaquettes alone [12] 



F = 



W 



lxl 



w 



lxlt 



Traceless 



(6) 



where W lxl is the clover-sum of four lxl Wilson loops lying in the //, v 
plane. For the three-loop improved topological charge operator, we employ 
an 0(a 4 ) improved definition of in which the standard clover-sum of four 
lxl Wilson loops lying in the /x, v plane is combined with 2x2 and 3x3 
Wilson loop clovers. Bilson-Thompson et al. [10] find 



-W 



lxl 



20w 4 



W 2X2 + 



90uf 



3x3 



h.c. 



(7) 



Traceless 



where W nxn is the clover-sum of four n x n Wilson loops and where F^ v is 
made traceless by subtracting 1/3 of the trace from each diagonal element 
of the 3x3 colour matrix. This definition reproduces the continuum limit 
with C(a 6 ) errors. On approximately self-dual configurations, this operator 
produces integer topological charge to better than 4 parts in 10 4 . 

We use periodic boundary conditions in space and anti-periodic in the time 
direction for fermions. Each of the hot configurations is also cooled for 12 
sweeps using a Cabibbo Marinari based algorithm in which the three diagonal 
SU(2) subgroups of SU(3) are looped over twice [12] using a three-loop C(a 4 ) 
improved action [10,11]. Twelve sweeps of cooling is just enough to see clear 
structure on the topological charge and action densities, but it still preserves 
much of the original topology. To see correlations for some eigenmode den- 
sities we also need to use the topological charge density of 5-sweep cooled 
configurations. These topological charge densities are much rougher, but are 



4 



closer to the original configuration, with topological objects moving less than 
in the 12-sweep cooled configurations. 

The visualisations of the eigenmode probability densities are treated consis- 
tently in order to allow direct comparison. This means they are all normalised 
with respect to the maximum value of each eigenmode. In this way, an isosur- 
face at half the peak height will reflect the size of the object. 



3 Smooth instanton background 



We generate a single instanton background on an 8 3 x 16 and on a 16 4 lattice 
by performing the path integration of 



M x ) 



X 



x 2 + p 2 \g 



(d,S)S-\ 



(8) 



with 



S = 



x 4 



(9) 



where + is an instanton and — an anti-instanton, to create the link variable. 
We find in the regular gauge 



u u g ( x ) = ex P io-nix) ■ (x; p) 



(10) 



<i>»(x;p) 



tan 



P 2 + ^i/jtfii-Ev X u ^ 2 



yV + Ev^u - X U ) 2 " P 2 + Ev(Xu ' X V ) 2 + {X V - X V ) 



ai(x) = (-X4 + x 4 , x 3 - x 3 , -x 2 + x 2 ), 
a 2 {x) = {—xz + X 3 , —X4 + X 4 , Xi - Xi), 
a 3 (x) = (x 2 - x 2 , -x 1 + xi, -x 4 + x 4 ), 
a 4 (x) = (xi -x 1 ,x 2 - x 2 , x 3 - x 3 ). 

In the singular gauge we find 6 

U^(x) = expK(x) • a(^(x; 0) - M (x; p)) 



(12) 



(13) 



6 Note that the analogous result of [3] inverts the roles of instantons and anti- 
instantons 
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bi(x) 
b 2 (x) 
b 3 (x) 
6 4 (x) 



(x 4 - x 4 , x 3 - x 3 , -x 2 + x 2 ), 
(-x 3 + x 3 , x A - £4, Xi - Xx), 
(x 2 - X 2) —Xi + X\, X4 - x 4 ), 



(-Xi + Xi, -x 2 + x 2 , -x 3 + x 3 ). 



(14) 



The singular gauge instanton is clearly recognisable on the volume rendered 
action density plot as seen in Fig. 1 (a). The outer surface shown is half the 
peak height. 

The results seen on the smaller exploratory lattice are also found on the 16 4 
lattice. The instanton on the latter lattice is cooled for two sweeps to minimise 
boundary effects. Eigenmodes of this configuration are calculated for 0.12 < 
k < 0.27 with an increment of 0.01 between values. The first four eigenvalues 
of the spectrum are shown in Fig. 2. We then evaluate the localisation of p(x) 
for the first three low-lying eigenmodes for each k. This is done by plotting 
p{x) as seen in Fig. 1 for selected eigenmodes. 

For eigenmodes with k <C k c — 0.125 we find a rather uniform distribution of 
p(x) whereas eigenmodes for all k c < k < 0.27 are strongly localised. For the 
lowest eigenmode this localisation strongly corresponds to the localisation of 
the instanton in the action density plot as seen by comparing Fig. 1 (a) and 
Fig. 1 (b), (c). We therefore say the eigenmode displays strong correlation 
with the instanton. The size and shape of this correlation varies from broad 
with some wall like structures for k < k c to very small for k c < k < 0.19 
and is getting broader again for larger values of k. For higher eigenmodes this 
localization gets broader and less correlated with the instanton. Some higher 
eigenmodes show no correlation with the instanton but just wall like structure 
as seen in Fig. 1 (f). Others show prolate-like correlations as seen in Fig. 1 (d). 
Some eigenmodes also show a wall-like structure and a prolate-like correlation 
as seen in Fig. 1 (e). It is useful to note that plane wave behaviour would 
display uniform behaviour in the density plot. 



4 Monte-Carlo generated SU(3) gauge field background 

After seeing that the low-lying eigenmodes of the Wilson-Dirac operator are 
strongly correlated to the instanton on a one-instanton configuration, we next 
carry out investigations on standard Monte-Carlo generated SU(3) background 
fields. We investigate three hot and three cooled configurations on both 8 3 x 16 
(f3 = 4.38) and 16 3 x 32 (/3 = 4.60) lattices. Three-loop improved cooling 
[10,11] is used to cool the configurations for 12 sweeps. This is just enough 
to get clear structure in the action and topological charge. As a reference for 
comparison of the eigenmode density we use the action and topological charge 
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Fig. 1. (a) Action density of a single instanton configuration on a f6 4 lattice, (b) 
First eigenmode for the single instanton configuration for k = 0.19. Note the very 
strong correlation with the instanton on the action density, (c) First eigenmode 
for the single instanton configuration for k = 0.25. Again note the strong correla- 
tion between the eigenmode and the action density, (d) Second eigenmode for the 
single instanton configuration for k = 0.25. The localization has a prolate shape 
compared to the spherical instanton. (e) Second eigenmode for the single instanton 
configuration for k = 0.19. The localization has a wall like shape with a prolate cor- 
relation to the instanton. (f) Third eigenmode for the single instanton configuration 
for k = 0.19. The eigenmode has a wall like shape and is not correlated with the 
instanton. For all figures what is shown is a volume rendering of the corresponding 
density. The outer isosurface is half the peak density. 
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Fig. 2. Smallest four eigenvalues of the eigenvalue spectrum of a smooth single 
instanton configuration on a 2-sweep cooled 16 4 lattice. For the open triangles we 
have strong spherical correlation between the eigenmode and action density and 
behaviour as described in the text. For the other symbols we have correlation with 
prolate- like shapes and wall-like structures. The shapes are shown in Fig. 1. For 
k = 0.20, 0.21 and 0.22 we expect the fifth eigenvalue to lie degenerate with the 
open triangles. 



density of the cooled configuration. We do this because those densities are too 
rough for the hot configuration and no accurate comparison would be possible. 
As there is more structure observed on the topological charge density plot, we 
use this as our preferred reference. See Fig. 4 (a) for a typical topological 
charge density plot of a 16 3 x 32 configuration. However the action density 
plot is also useful for guidance. 

We say an eigenmode is correlated to a topological object if p(x) has a peak 
within one lattice site of the peak topological charge density. We calculate the 
first four eigenmodes for values of 0.13 < n < 0.25 for cooled configurations 
and 0.15 < k < 0.29 for hot configurations in steps of 0.01. This range is from 
approximately 0.02 smaller than k c to a region where doublers appear in the 
overlap formalism. We find that the behaviour described in the following is 
general for all configurations: Each of the eigenmodes is localised. This was 
already observed for the lowest eigenmode [5]. For k < k c this localization 
weakens and the density, p(x), broadens quickly. An exponential decay of the 
density p(x), as previously observed [2,5], seems likely to occur. All low-lying 
eigenmodes for k > correlated to topological objects. 
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Fig. 3. Smallest four eigenvalues of the eigenvalue spectrum of a hot configuration 
showing correlations of eigenmodes to topological objects. Each eigenmode has a 
correlation to one topological object and each different symbol indicates a different 
topological object. Modes for k = 0.15 < k c are very broad and weakly correlated 
to several objects. 



We are able to track a correlation along the modes for all k > k c and for 
one step smaller than k c . Thus we can label a mode by its correlation. Fig. 
3 shows eigenmodes of a hot configuration, where the symbol used denotes 
which topological object the localised low-lying eigenmode is correlated with. 
Careful inspection of Fig. 3 reveals the presence of an eigenmode, correlated to 
one topological object, but with a spectral flow containing two zero crossings. 
The size of the correlated object, obtained in the following section, suggests 
that this object is a lattice artifact. 

Around k c and towards the upper end of the analysed spectrum, localisations 
are weak. However, the correlation of positions is strong. At these k we tend 
to get correlations to more than one topological object. The localisations get 
stronger and sharper for increasing values of k until a maximal localization is 
reached for a value k we will define to be K max . For k > /t max the localisations 
get weaker again. Fig. 4 shows a visualisation of this behaviour. Additional 
figures can be found elsewhere [13]. We take a closer look at this behaviour in 
the next section. 



For the lower eigenmodes, which are separated by a gap from the higher eigen- 
modes, we find correlation to one topological object per eigenmode. For higher 
eigenmodes , which are closer together with some degenerate modes, we find 
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Fig. 4. (a) Topological charge density of a 12-sweep cooled 16 3 x 32 configuration at 
[3 = 4.60. (b) Action density of the same configuration. There is clearly less struc- 
ture for the action density than the topological charge density, (c) The low-lying 
eigenmode density at k = 0.14 of the same configuration with correlation to the 
object seen in the action density and the topological charge density, (d) Eigen- 
mode density at k = 0.16 showing correlation to the same topological object, (e) 
Eigenmode density at k = 0.18. The localization is maximal for this value of k. (f) 
Eigenmode density at k = 0.19. The localization is becoming larger again. 
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correlations to more than one topological object. In general we find that the 
closer the eigenvalues and the weaker the actual localisations, the more likely 
it is to get correlation with more than one topological object in the corre- 
sponding eigenmodes. 

Calculations of 20 eigenmodes for selected values of k on hot configurations 
show only little broadening in the localisations and persistent correlations with 
topology. This suggests that this behaviour will not change quickly and such 
correlations will persist for even higher eigenmodes. 

As mentioned above, we use the topological charge density of the 12-sweep 
cooled configuration for comparison. The position of eigenmodes on cooled 
fields agree perfectly well with the position of topological objects seen in this 
density. A typical distance between such positions being 0.05 lattice units. For 
hot configurations some of the correlations are slightly offset compared to the 
structure in the smoothed configurations whereas other structure in the eigen- 
modes has no corresponding structure in the smoothed topological charge. But 
in those cases a comparison with a less cooled configuration reveals correlations 
between eigenmodes and topological objects which are moved or destroyed by 
further cooling. This behaviour is expected and understood, as topological 
objects are known to move under cooling as instantons and anti-instantons 
attract each other and annihilate when they are close enough together. 

This clear correlation between p(x) of eigenmodes on hot configurations and 
topological objects suggests that it is possible to identify areas through the 
noise of a hot configuration with significant topology. We can "see through" 
the noise by using eigenmodes of the hermitian Wilson-Dirac operator. 

The spectral plot of the first four eigenvalues of a 12-sweep cooled configuration 
looks different compared to the same plot on a hot configuration. Comparing 
Fig. 3 and Fig. 5 we notice the rhomboid shape with an area without any 
eigenvalues in the spectral plot of the cooled configuration. This is expected 
for such smoothed configurations [14]. Eigenmodes on the right-hand side of 
this rhomboid behave different than eigenmodes on the left-hand side of the 
rhomboid. The localization of the eigenmodes on the left-hand side is approx- 
imately the same as the localisations in the hot configuration. But the eigen- 
modes on the right-hand side are much weaker in localization. Those weaker 
localised eigenmodes are related to very high eigenmodes in the hot configura- 
tion and the process of cooling brings them into an area where we can observe 
them as low eigenmodes. These eigenmodes must be very high eigenmodes 
in a hot configuration for they are much weaker localised as localisations we 
have seen for up to the 20th eigenmode in a hot configuration. Although the 
behaviour in localization strength for eigenmodes of cooled configurations is 
different from that of eigenmodes of hot configurations the correlations with 
topological objects still exist for all eigenmodes. We will take a closer look at 
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Fig. 5. First four eigenvalues of the spectrum of a 12-sweep cooled configuration. 
Note the rhomboid shape with the area without eigenvalues. The modes to the right 
of the maximum, the filled symbols, are only weakly localised. The open symbols 
are strongly localised. 



localization strength and quantify it in the next section. 



5 Quantitative Results 



As described earlier in this paper, the localization of the eigenmodes change 
shape and size with changing of k. In order to quantify this behaviour we 
have two methods. The first one is fairly simple. We have already seen that 
most eigenmodes are localized at one topological object. As the eigenmodes 
are all normalized, the maximum, or peak, value of the eigenmode density is 
an indicator of how strong this localization is. Fig. 6 shows the plot of such 
peak values for four eigenmodes each of three hot 16 3 x 32 configurations. 
This plot shows a smooth behaviour with a maximum for /t max ~ 0.23 for 
the lowest eigenmodes. This suggests that for hot SU(3) configurations the 
strongest localisation occurs in that k region. 

Fig. 7 shows the same plot for four eigenmode densities each for three 12-sweep 
cooled 16 3 x 32 configurations. This plot shows a different behaviour. Again 
there is a clear peak, but the values drop to almost zero immediately after 
the peak. We expect this to happen after our qualitative observations in the 
previous section. These low maximum values correspond to weakly localized 
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with respect to k 



0.07 
0.06 
0.05 f 

0) 

a 0.04 h 
> 

ft 

0.02 
0.01 - 
0.00 



-X L 



X 



X 



X 





X 


X 


X 




X 


X 


X 


X 




X 


X 




X 




X 


X 





X 



J X- 



x 



X 



X 



X 



X 



-x^ — * — H — x- — x- 



0.12 0.14 0.16 0.18 0.20 0.22 0.24 0.26 

Fig. 7. Peak value of four eigenmode densities each from three 12-sweep cooled 
16 3 x 32 configurations with respect to k 
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modes, which are on the right-hand side of the spectrum shown in Fig. 5. In 
this case we can not determine K max , for it is not clear how the localisations are 
going to develop if we follow the stronger localised eigenmodes up to higher 
values of n. The one eigenmode we could follow suggests that 0.20. 
This reduction of /t max is about the same size as the reduction of k c for going 
from the hot to the 12-sweep cooled configuration. 

To get more information about the shape of the eigenmodes, we fit the in- 
stanton model, Eq. (4), to the eigenmodes. This model gives us very good fits 
as 



where Pm{x) is the 6 parameter fit of 3 4 points of p(x) centered about the peak 
of p(x). This value is about 1000 times smaller than the peak value and 10~ 5 
is the worst case with most of the fits of order 10~ 7 to 10 -9 . The fit parameter 
p is then a good measurement for the size of the localization. 

In Fig. 8 we plot p(n) for four eigenmodes calculated on all hot configurations 
and in Fig. 9 for four eigenmodes calculated on all 12-sweep cooled configura- 
tions. We see a behaviour which corresponds to the behaviour described above 
for the peak values of the eigenmodes. p(k) for the hot configurations shows 
a smooth behaviour with a minimum around « max = 0.26. This means the 
eigenmodes are maximally localised for this value of k and are less localized 
for both ends of the spectrum. The K max found this way varies slightly from 
K max found by just taking the peak values of the eigenmode densities. It is 
about 0.03 larger at 0.26 for the lowest eigenmode. As we mentioned earlier 
the eigenmodes for low and high values of k are localized on more than one 
topological object. Therefore Fig. 8 reports more local maxima than modes 
at these k. However in the range 0.18 < k < 0.28 only one local maximum is 
found per mode. 

We have established that low-lying eigenmodes are correlated to a single topo- 
logical object when eigenmodes are non degenerate. It is also established that 
an instanton gives rise to a zero crossing in the spectral flow with the sign of 
the slope equal to the sign of the topological charge [3]. However, as mentioned 
above, Fig. 3 reveals an eigenmode correlated to one topological object, but 
with two zero crossings. The size of the eigenmode varies from p = 1.2 to 
p = 1.0 as k varies from 0.21 to 0.27. The size of the correlated object on the 
topological charge density, which can only be seen on the 5-sweep-cooled con- 
figuration, is p = 1.1. We have been able to reproduce similar spectral flows on 
single instanton configurations, which have been cooled with the Wilson ac- 
tion to the point where the topological charge is Q ~ 0.4, well below 1. Hence 
the double zero crossing of the spectral flow in Fig. 3 suggests the presence 
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Fig. 8. p of the fitted instanton model of four eigenmodes from three hot 16 3 x 32 
configurations with respect to k. Note that there are several correlations in one 
eigenmode at both ends of the spectrum. 

of a lattice artifact. Ideally an improved fermion action should act to remove 
the zero crossings associated with this artifact. 

The behaviour of p(n) of the cooled configurations shows a jump at k m 0.19, 
where the weakly localised modes set in as described above. Some of these 
weakly localised modes are not sufficiently localised to allow a fit to the instan- 
ton model. A fit would result in values for p > 10, which are not reasonable for 
the instanton model and therefore neglected. For those eigenmodes where we 
can do a fit, we get large values of p compared to the strongly localised eigen- 
modes. We find that p is about 2 for strongly localised eigenmodes and about 
6 for weakly localised eigenmodes. For k < 0.19, p(/t) decreases smoothly, 
but for k > 0.19 the values of p(k) are higher than at the lower end of the 
spectrum and do not show a smooth behaviour. Again it is hard to extract 
K max , but it seems that if we could follow the stronger localised modes further 
it would be around 0.22. This is about 0.02 larger than the « max , which is 
extracted using the peak values. 

The fitting of the models to the eigenmode densities, using Eq. (4), as well 
as to the topological charge densities, using Eq. (3), allows us to compare the 
sizes for the eigenmode localisations with the sizes of the actual topological 
objects. In order to do this, we find the topological object which is located 
closest to the position of the eigenmode. We can fit the instanton model only to 
topological charge densities of cooled configurations, as only those are smooth 
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Fig. 9. p of the fitted instanton model of four eigenmodes from three 12-sweep cooled 
16 3 x 32 configurations with respect to k. Note the jump at k « 0.19. 

enough. Therefore we only compare the sizes of eigenmodes of cooled config- 
urations with sizes of actual topological objects. Due to the different localiza- 
tion strength for eigenmodes on the left-hand and on the right-hand side of 
the rhomboid spectrum of a cooled configuration, we look at those eigenmodes 
separately. 



We find that the fitted positions for the strongly localised eigenmodes agree 
very well. With the fitted positions of the correlated topological objects lie 
within a fifth of a lattice spacing, and most of the times even better. The 
fitted positions of the weakly localised modes to the right of the spectrum 
agree only within one lattice spacing with their correlated topological objects. 
Strongly localised eigenmodes are correlated to smaller objects in the topolog- 
ical density. Such eigenmodes on the left-hand side of the eigenvalue spectrum, 
except those for k < k c , have a size between p pa 1.5 and p pa 3.5. The corre- 
lated topological objects have a size between p pa 2 and p pa 3. The eigenmodes 
are larger than their correlated topological objects for smaller values of k, but 
as they shrink with growing k they get smaller than their correlated topolog- 
ical objects. All followed modes reach the size of the correlated topological 
object for 0.155 < k < 0.175. Fig. 10 shows an example of this behaviour. 

The weakly localised eigenmodes on the right-hand side of the eigenvalue spec- 
trum are larger than the strongly localised eigenmodes on the left-hand side of 
the spectrum. They turn out to have a size between p pa 4 and p pa 6. But the 
topological charge density correlated with those eigenmodes are themselves 
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Fig. 10. p of the fitted instanton model of one mode of a 12-sweep cooled 16 3 x 32 
configuration with respect to k. The x denote the size of the eigenmode and the + 
denote the size of the correlated topological object. That the two graphs cross is a 
general result and all crossings are found for 0.155 < k < 0.175. 

bigger than the topological charge density correlated with the stronger lo- 
calised modes. The sizes of the topological charge density objects lie between 
p PS 3 and p m 4. The weaker localised eigenmodes are always bigger than the 
correlated topological objects with p mo dc ~ PtopQ + 2. 

For k — k c the size of the eigenmodes is between p = 3 and 5 with a relation 
to the size of the correlated topological objects of p mo dc ~ PtopQ + 1- Sizes of 
strongly localised eigenmodes of hot configurations are of a comparable, but 
slightly smaller size than that on cooled configurations. 

In order to understand why modes on the right-hand side of the rhomboid are 
just weakly localised we calculate spectra for one configuration with different 
amounts of cooling. We find that the value of k where the jump in the locali- 
sation size occurs becomes smaller with cooling. Compare Fig. 9 and Fig. 11. 

It is known [4] that zero crossings for larger k correspond to smaller topolog- 
ical objects. As we see localisation sizes of calculated eigenmodes shrink with 
growing k as well, we generally think of larger values of k being associated 
with smaller objects. Cooling removes smaller objects first, therefore eigen- 
modes for large n become "unassociated" with small objects. The eigenmodes 
show a behaviour of weak localisation which can be seen in Fig. 9 and Fig. 11 
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Fig. 11. p of the fitted instanton model of four eigenmodes for one 4-sweep cooled 
16 3 x 32 configuration with respect to k. Note the jump occurring at higher k 
compared to Fig. 9. 

after the jump. Further cooling removes larger objects and therefore the value 
of k, where the change of behaviour sets in, becomes smaller. 



6 Conclusion 



Not only zero modes, but all low-lying eigenmodes of the hermitian Wilson- 
Dirac operator Dw are strongly correlated to topological objects for k > k c . 
These objects can be instantons for which S/So = Q or topological fluctuations 
with S/So > Q. One eigenmode is correlated to at least one topological object 
with correlations to more than one object becoming more likely as eigenvalues 
become degenerate and eigenmodes become broader in size. For k < k c the 
correlations broaden very quickly and are lost for values smaller than about 
(k c — 0.02). For k > k c the correlations become sharper until n = K max for 
which the correlations are strongest. For k > /t max the correlations broaden 
again. 

Eigenmodes of 12-sweep cooled configurations show a different behaviour de- 
pending on whether they belong to eigenvalues on the left or right-hand side 
of the rhomboid-shaped eigenvalue spectrum of a cooled configuration as seen 
in Fig. 5. Eigenmodes belonging to the left-hand side of the eigenvalue spec- 
trum are strongly localised and show the same behaviour as eigenmodes of hot 
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configurations. Eigenmodes belonging to the right-hand side of the eigenvalue 
spectrum are very weakly localised, but are still correlated to topology. This 
suggests that those eigenmodes correspond to very high eigenmodes in a hot 
configuration, which are lowered by cooling. 

The value of k where the weakly localised eigenmodes set in becomes smaller 
with cooling. This supports the idea of high values of k corresponding to lo- 
calisations on small topological objects. Small topological objects are removed 
first under improved cooling thus eigenmodes for high values of k are the first 
ones to loose the strong localisation. 

When an instanton model is fitted to the eigenmode density, using Eq. (4), 
and to the topological charge density, using Eq. (3), strongly localised eigen- 
modes have about the same size for p as the correlated topological objects. 
For n max eigenmodes are slightly smaller than the correlated topological and 
for k smaller than — 0.02 eigenmodes are slightly larger. 

On a single instanton configuration the correlation to the instanton persists 
strongly only for the lowest eigenmode and is then gradually lost for higher 
eigenmodes. On SU(3) background fields we see correlation for the 20 lowest 
eigenmodes. There is only little broadening which suggests that the correlation 
will persist for eigenmodes higher than 20. 

For hot configurations these correlations allow us to "see through" the noise 
to underlying topological objects. This enables one to track the movement of 
these objects as a function of cooling. 

Future work will examine the manner in which the eigenvalue spectrum is 
modified under improved fermion actions. In particular we plan to examine the 
effects of using APE-smeared fat-links in the irrelevant operators of fermion 
actions. 
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